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By analyzing chaotic states of the one-dimensional 
Kuramoto-Sivashinsky equation for system sizes L in the 
range 79 < L < 93, we show that the Lyapunov fractal di- 
mension D scales microextensively, increasing linearly with L 
even for increments AL that are small compared to the av- 
erage cell size of 9 and to various correlation lengths. This 
suggests that a spatially homogeneous chaotic system does 
not have to increase its size by some characteristic amount 
to increase its dynamical complexity, nor is the increase in 
dimension related to the increase in the number of linearly 
unstable modes. 
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An important phenomenon associated with sustained 
nonequilibrium systems is spatiotemporal chaos, a 
chaotic dynamical state that is spatially disordered . 
An open question is how best to characterize spatiotem- 
poral chaos so that theory can be quantitatively com- 
pared with experiment and experiment with simulation. 
Presently, there is no fundamental theory of nonequilib- 
rium systems to indicate the appropriate quantities to 
measure and so researchers have borrowed ideas from 
condensed matter physics, fluid dynamics, nonlinear dy- 
namics, and statistics. Commonly used ways to charac- 
terize spatiotemporal chaos include critical exponents || , 
the two-point correlation time T2 0), the largest Lya- 
punov exponent Ai (2), the Lyapunov fractal dimen- 
sion D, the two-point correlation length £2 H, the di- 
mension length £s HQ and other lengths (|,|7]-|9|] . How- 
ever, calculations have shown that these quantities do not 
always lead to the same conclusions, e.g., there are sys- 
tems for which the length £2 diverges while the length £5 
remains finite as some parameter is varied 0. Further 
research is therefore needed to understand the particular 
features of spatiotemporal chaos that are measured by 
any one of these quantities and how these quantities are 
related to one another. 

In the following, we report results that provide new in- 
sights about how the dynamical complexity of a nonequi- 
librium system depends on the volume of the system, 
and about the interpretation of the dimension length £5. 
In 1982, Ruelle conjectured [10j, and numerical calcula- 



tions later confirmed 1 1 1 



that the dimension D 
of a sufficiently large spatially homogeneous chaotic sys- 
tem should increase extensively, i.e., linearly with its vol- 



ume V. Using an argument similar to that used by Lan- 
dau and Lifshitz to explain the extensivity of additive 
quantities in thermodynamics |l5| ], this extensivity of D 
can be understood heuristically as a consequence of spa- 
tiotemporal disorder. If two subsystems of a spatiotem- 
poral chaotic system are sufficiently far apart, their cou- 
pling is weak because of the disorder and so their dynam- 
ics contribute independently and additively to the overall 
fractal dimension. 

This picture of weakly interacting subsystems raises 
the question of how precisely does the fractal dimen- 
sion D increase with increasing system volume in the 
extensive regime. One possibility is that the curve D(V) 
may be linear only on average and has a staircase-like 
structure, with the steps corresponding to new degrees 
of freedom that appear once the system volume has in- 
creased sufficiently to include a new subsystem. The 
widths Ay of the steps would then define a length 
scale (AV) 1 ^ (where d is the spatial dimensionality of 
the system) that would be interesting to compare with 
the lengths mentioned above (£2, £5, etc.). Possible step- 
like features in the D(V) curve might also be associated 
with the appearance of new linearly unstable modes of 
the uniform state, since the number of such modes typ- 
ically increases linearly on average with increasing vol- 
ume p6| . Another possibility is that the curve D(V) is 
extensive only on average but its deviation from linearity 
is too irregular to characterize by a single length scale. A 
fourth possibility is that there are no length scales asso- 
ciated with how D increases with V and the curve D(V) 
is exactly linear for arbitrarily small increases in V, a sit- 
uation that one could call microextensive chaos. In this 
case, it would be interesting to understand how the geo- 
metric structure of the chaotic attractor in phase space 
changes with V so as to produce such an exact linear 
behavior. 

In this paper, we numerically integrate the one- 
dimensional Kuramoto-Sivashinsky (KS) equation — a 
widely studied continuum model of spatiotemporal 
chaos H — to investigate how the Lyapunov fractal di- 
mension D(L) of a homogeneous chaotic system varies 
with the system size L for increments AL that are small 
compared to the lengths mentioned above (£2, £,s, etc). 
With one exception Q , all prior numerical studies used 
increments AL that were large compared to these lengths 
and the detailed form of D(L) was not determined. We 
show below that, in fact, the Lyapunov fractal dimen- 
sion D increases linearly with L even for system incre- 
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ments AL that are tiny compared to the average cell size 
and to various correlation lengths. The spatiotemporal 
dynamics of the one-dimensional KS equation therefore 
provides an example of microextensive chaos. We con- 
jecture that this will be a general property of chaotic 
homogeneous nonequilibria media. 

Our calculations of D versus L yield an additional in- 
sight, namely that the onset of extensivity in D is not 
sharp but occurs only asymptotically with increasing L, 
after a sequence of alternating windows of stationary, pe- 
riodic, intermittent, and chaotic dynamics. (Such alter- 
nating windows have been noted before |f7j but have 
not been studied with such fine resolution in L as we do 
here.) These results suggest the possibility that windows 
of non-chaotic behavior may persist to arbitrarily large 
values of L but become too narrow to be detected. If 
true, then the dimension D(L) may not be a continuous 
curve and extensive behavior occurs only between the 
narrow windows of non-chaotic dynamics. 

Our results were obtained by numerical integrations of 
the onc-dimcnsional Kuramoto-Sivashinsky equation in 
the form 

d t u(t,x) — —d x u — d*u — ud x u, i 6 [0,1], (1) 

on an interval of length L, with rigid boundary condi- 
tions u — d x u — at x — and at x = L, (Fig. [| shows 
a chaotic and periodic state for L = 50 and 54 respec- 
tively.) The spatial derivatives were approximated by 
second-order-accurate finite differences on a uniform spa- 
tial mesh, and a standard operator-splitting method was 
used for the time integration [ p~8| . For given initial con- 
ditions and interval length L, we used the Kaplan- Yorke 
formula pl| to calculate the Lyapunov fractal dimen- 
sion D(L) in terms of all of the positive and some of the 
negative Lyapunov exponents A^. These exponents were 
obtained using a standard algorithm |l9tl in which many 
copies of the linearized KS equation were integrated, each 
with their own initial condition. 

The demonstration of microextensive scaling by the 
above numerical methods was delicate since the Lya- 
punov exponents A.; and so D converge noisily and 
slowly pej ] toward their infinite time limits. As the in- 
crement AL in system size became smaller, the corre- 
sponding increment in dimension AD was more difficult 
to determine since AD became comparable to the fluc- 
tuations in the dimension curve D(T) as a function of 
integration time. The exponents A.; and dimension D(L) 
were also sensitive to the values of the spatial resolu- 
tion Ax and temporal resolution At, to the renormaliza- 
tion time Tnorm for the Lyapunov vectors, and to the total 
integration time. For nearly all runs reported below, we 
used values of Ax = 0.167, At = 0.025, and T norm = 10 
and confirmed the correctness of the corresponding re- 
sults by comparing the values with spatial and temporal 
resolutions up to four times larger and for integration 
times as long as 10 6 time units. 



We now turn to our results. Our starting point was the 
pioneering calculation of Manneville fil| ] , who used nu- 
merical integrations of Eq. (|l]) with rigid boundary con- 
ditions to demonstrate for the first time that the fractal 
dimension D scaled extensively with the system size L. 
For L > 50, he found that D = 0.230L - 2.70, which 
implies a dimension length Q of = (dD / dL)~ l = 
1/0.230 w4.4@, This length is somewhat smaller than 
the average cellular size A = 27r/<7 max = 2\/2tt w 8.8 cor- 
responding to the fastest growing linear mode g max = 
l/-\/2. Based on these results, we chose to calculate the 
fractal dimension D{L) over the range 50 < L < 100 in 
constant increments AL — 0.5 that were much smaller 
than these lengths. In contrast, the smallest increment 
used by Manneville was AL — 50 for which the fractal 
dimension changes by about 12. 

Manneville's linear dependence of D on L suggested 
that for L > 50, only spatiotemporal chaos exists. In 
contrast, we find that there is a complicated sequence of 
different dynamical states over the range 50 < L < 75 
and then only chaotic states for 75 < L < 93 |^3| . Fig. || 
summarizes our results for the range 50 < L < 75 by 
plotting the period of each state as a function of L. We 
observe four kinds of states: fixed points, time-periodic 
states, chaos, and intermittent states in which one kind 
of time dependence alternates irregularly with a different 
kind of time dependence. In most cases, these different 
categories were easily identifiable to the eye by looking 
at time series. To combine all the results on a single plot, 
chaotic states were arbitrarily assigned a period of -200, 
intermittent states a period of -100, fixed points a period 
of 0, and periodic states a direct estimate of their period 
based on repeating features of the time series. 

There are two interesting features of the dynamical 
states of Fig. |^ in addition to the unexpected occurrence 
of many windows of alternating dynamics for this range 
of L. First, we found that for a given system size L, 
numerical integrations using up to seven different ran- 
dom initial conditions (each consisting of uniformly dis- 
tributed numbers in the interval [—0.1,0.1]) led to only 
one state. Thus empirically there seems to be only one 
basin of attraction for each system size and we do not 
expect hysteresis in the range 50 < L < 75. Second, we 
found rather remarkably that the fractal dimension D of 
each chaotic state in Fig. ^ lay on Manneville's exten- 
sive curve D(L) = 0.230L - 2.70 with D > L>(50) = 8.8 
(A least-squares fit of our chaotic states gave the almost 
identical curve 0.227L — 2.85). Thus the states jump 
abruptly from low-dimensional D = 1 periodic states to 
high-dimensional chaotic states that are scaling exten- 
sively with the system size. We did not try to character- 
ize the intermittent states, e.g., by their fractal dimension 
or by the scaling properties of the fractional duration of 
a particular phase pl[ . 

Over the range 78 < L < 93, only chaotic states were 
observed. Fig. S shows that the corresponding values 
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of the Lyapunov fractal dimension D lie on a straight 
line that has the same slope (to two digits) and inter- 
cept as that found by Manneville over the much larger 
range 50 < L < 400. Thus the fractal dimension shows 
microextensive scaling: a linear dependence on L for sys- 
tem increments AL — 0.5 that are much smaller than any 
characteristic length scale such as the average cell size or 
various correlation lengths. Given the similar result ob- 
tained by Xi et al (l| for a different mathematical model 
of spatiotemporal chaos, we conjecture that microexten- 
sive scaling will be a general feature of spatiotemporal 
chaos in sufficiently large, approximately homogeneous 
nonequilibrium systems. 

The linear behavior of Fig. [| rules out a simple relation 
between the Lyapunov fractal dimension and the number 
of linearly unstable modes. For both periodic and rigid 
boundary conditions, the wave numbers are quantized 
in units of 2tt/L and so the number of unstable modes 
increases linearly as L/tt « 0.32L, but in discrete jumps 
when L changes by about it. Over the range 78 < L < 93, 
we would expect (93 — 78) /it ss 4 new linearly unstable 
modes to appear but there are not correspondingly four 
step-like features in the figure. 

Fig. ^ shows how the first sixteen exponents Xi vary 
with system size L over the same range of system sizes 
as Fig. ||. The largest exponent Ai is approximately inde- 
pendent of L and so is an intensive quantity; an increase 
of system size therefore does not change the forecasting 
time 1/Ai over which the future behavior of the field u 
becomes unpredictable, despite the fact that the fractal 
dimension is becoming correspondingly larger. The other 
exponents increase roughly linearly with increasing L, 
with the negative exponents having the largest slope and 
showing substantial deviations from a strictly linear de- 
pendence on L. Comparing Fig. |]with the Kaplan- Yorke 
formula for the dimension 
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where K is the largest integer such that 5D i=1 \ > 0, we 
see that the linear behavior of D with L in Fig. ^ is a 
subtle consequence of how the Ai vary with L. 

In conclusion, we have demonstrated the occurrence 
of microextensive scaling of the Lyapunov fractal dimen- 
sion with system size for the one-dimensional Kuramoto- 
Sivashinsky equation with rigid boundary conditions. 
This suggests that a spatially extended nonequilibrium 
dynamical system does not have to increase its volume 
by some minimal amount for the fractal dimension to in- 
crease. Correspondingly, the dimension length £5 does 
not have some direct physical meaning as the character- 
istic size of a dynamical subsystem, it is simply deter- 
mined by the linear growth of D with L. Our calcula- 
tions suggest two questions for further exploration. One 
is whether there is a cutoff system size L c above which 



only chaotic solutions are found for the KS eq. Second 
is to understand mathematically how the geometry of a 
strange attractor changes with system size L such that 
that D varies exactly linearly with L. 
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FIG. 1. Space-time evolution of the field u(t, x) for 
two states of the Kuramoto-Sivashinsky equation Eq. (|l]) 
with rigid boundary conditions. The space-time resolution 
was At = 0.025 and Ax — 0.166 and the peak-to-peak am- 
plitude is about 4. (a): Chaotic state for L — 50. Spa- 
tial curves are plotted every AT = 1 time units starting at 
time t — 50,000. (b): A periodic state for L = 54 with pe- 
riod r = 127.6. Spatial curves are plotted every AT = 5 time 
units starting at time t = 80, 000. 
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FIG. 2. Periods of numerical solutions to Eq. QlJ) ver- 
sus the system length L. Each integration was started from 
small-amplitude random initial condition and then integrated 
500,000 time units. There is a complex sequence of windows 
corresponding to chaotic, constant, intermittent, and periodic 
dynamics. Chaotic and intermittent solutions have been as- 
signed an arbitrary period of -200 and -100 respectively so 
that all the data could be compared on one plot. 
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FIG. 3. The Lyapunov fractal dimension D of chaotic so- 
lutions to Eq. ([l]) versus system size L for 79 < L < 93. The 
dimension values accurately fall on a straight line, demon- 
strating the occurence of microextensive scaling. The straight 
line D — 0.227L — 2.85 was obtained by a least-squares 
fit to the points and agrees well with Manneville's re- 
sult of D — 0.230-L — 2.70 over the much larger range 
50 < L < 400. The error bar for each point corresponds to 
a relative error of at most 0.05% in D. The error bar was 
determined by the peak-to-peak fluctuations of D versus in- 
tegration time T. 
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FIG. 4. The sixteen largest Lyapunov expo- 

nents Ai > ■ ■ ■ > Ai6 versus system length L for the same 
window of chaotic solutions discussed in Fig. ^. For each value 
of L, the Ai were estimated by estimating the asymptotic slope 
of the Ai(t) versus t curve for 20, 000 < t < 200, 000. All Lya- 
punov exponents increase with L in this regime except Ai, 
which is approximately independent of L and so is an inten- 
sive quantity. The lines through the points for a given Ai are 
a guide to the eye and were determined by a least-squares fit. 
The exponents Ai increases approximately, but not exactly, 
as linear functions of L. 
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